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Abstract 

We obtain the supergravity dual of A/" = 1 gauge theory in 2+1 dimensions with a large 
number of unquenched massless flavors. The geometries found are obtained by solving the 
equations of motion of supergravity coupled to a suitable continuous distribution of flavor 
branes. The background obtained preserves two supersymmetries. We find that when N^. > 
2Nf the behavior of the solutions is compatible with having an asymptotically free dual 
gauge theory with dynamical quarks. On the contrary, when < 2Nf the theory develops 
a Landau pole in the UV. We also find a new family of (unfiavored) backgrounds generated 
by D5-branes that wrap a three-cycle of a cone with G2 holonomy. 
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1 Introduction 



The gauge/gravity correspondence [H [2] has provided us with a very powerful tool to ex- 
plore the dynamics of gauge theories at strong coupling. In its original formulation the 
correspondence is a duality between the AdSt^ x background of type IIB supergravity 
and J\f = 4:,d = 4 super Yang-Mills theory, in which all fields transform in the adjoint 
representation of the gauge group. Clearly, to extend this duality to systems closer to the 
particle physics phenomenology we should be able to add matter fields transforming in the 
fundamental representation of the gauge group. This is equivalent to adding open string 
degrees of freedom to the supergravity side of the correspondence. 

It was originally proposed in [3J that such an open string sector can be obtained by adding 
certain D-branes to the supergravity background. If the number Nf of such flavor branes is 
small compared with the number of colors, we can treat the flavor branes as probes in 
the background created by the color branes. This is the so-called quenched approximation, 
which corresponds, in the field theory side, to suppressing quark loops by factors 1/Ncin the 
't Hooft large expansion. The fluctuation modes of the probe D-brane in the supergravity 
background provide a holographic description of the flavor sector of the gauge theory and one 
can extract the corresponding meson spectrum by analyzing the normalizable fluctuations 
of the probe [1] (for a review and a list of references, see [5]). 

When the number Nf of flavors is of the same order as the number Nc of colors, the 
backreaction of the flavor branes on the metric can no longer be neglected. On the field theory 
side the inclusion of the backreaction is equivalent to considering the so-called Veneziano 
limit, in which Nc and Nf are large and their ratio Nf/Nc is fixed. In this limit quark loops 
are no longer suppressed. In the last few years there have been several attempts to construct 
supergravity duals of these unquenched systems, both for four-dimensional [6] and three- 
dimensional gauge theories [7], by using solutions of supergravity generated by localized 
intersections of branes. 



Recently, a different approach has been proposed in ref. [8j. Instead of solving the equa- 
tions of pure supergravity, the authors of |8] considered the full gravity plus (flavor) branes 
system. The action of such a system contains the Dirac-Born-Infeld action of the flavor 
branes, which governs their worldvolume dynamics and their coupling to the different su- 
pergravity fields. Notice that this is consistent with the fact that color branes undergo a 
geometric transition and are converted into fluxes, whereas, on the contrary, flavor branes are 
still present after the geometric transition. Thus, from a conceptual point of view, color and 
flavor branes are not equivalent and, therefore, should be treated differently. By considering 
a suitable continuous distribution of flavor branes the authors of [Hj were able to find a set of 
BPS equations and to solve them numerically (see [9] for a similar approach in the context 
of non-critical string theory). The resulting solution is the flavored backreacted version of 
the background found in [10] and proposed in [llj as the supergravity dual of A/" = 1 super 
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Yang-Mills theory in four dimensions. Further developments of this approach can be found 
in refs. [I2]-[II]- In this paper we will apply this circle of ideas to the case of A/" = 1 gauge 
theories in 2+1 dimensions. The corresponding unfiavored supergravity dual was found in 
ref. jlH], and it was interpreted as being generated by D5-branes wrapped on a three-cycle 
of a manifold of G2 holonomy in [19] (see also [20]-[23]). 

The low number of supersymmetries preserved by the solution of (just two real super- 
charges) is a nice feature and makes it appealing also from the perspective of its dual field 
theory. As pointed out in ref. [19j, this theory reduces in the IR to 2+1 dimensional M = 1 
supersymmetric U (Ai"c)Yang-Mills theory with a level k Chern-Simons interaction. Such the- 
ory coupled to an adjoint massive scalar field should arise on the domain walls separating the 
different vacua of pure = 1 super- Yang-Mills in 3+1 dimensions. For 7^ the theory has 
a mass gap, at least classically, with mass of order gYM\^\- This implies that for \k\ » 1, 
i.e. when we can trust the classical result, there are no Goldstone fermions and, therefore, 
supersymmetry is unbroken. Actually, the Witten index for such a theory was computed 
in [21], where it has been shown that for k > Nc/2 supersymmetry is unbroken, while it is 
broken for k < Nc/2. In the borderline case {k = Nc/2) there is just one supersymmetric 
vacuum. Being the supergravity solution of [18] supersymmetric and without parameters 
that could label different vacua, it is reasonable to expect that the dual field theory is the 
one describing the k = Nc/2 case. It was shown in ref. [19] that this is actually the case. 

We will start our analysis by generalizing the ansatz of [18] for the unfiavored solutions. 
This generalization will allow us to find a new class of solutions in which, in the UV, the 
metric becomes asymptotically the direct product of a G2 cone and a three-dimensional 
Minkowski space, while the dilaton becomes constant. This is in contrast to the background 
of [18], in which the dilaton grows linearly with the holographic coordinate. For this gen- 
eralized ansatz we will be able to find a system of first-order BPS equations which ensure 
that our solutions preserve two supersymmetries. We will perform a careful analysis of the 
regularity conditions to be imposed on the functions of our ansatz, which will allow us to fix 
some parameters of our solutions and to determine the appropriate initial conditions needed 
to solve the BPS differential equations. The new solutions, which are found numerically, can 
be naturally interpreted as non-near horizon versions of the one of [T8] . 

After completing the analysis of the unfiavored backgrounds, we will study the addition 
of flavor D5-branes. First of all, we will use kappa symmetry |25| to determine a continuous 
family of embeddings of probes that preserve all the supersymmetries of the background 
and which can be used as flavor branes for massless quarks. These embeddings have the 
topology of a cylinder and are very similar to the ones found in [26] (and used in [8j) in the 
case of the supergravity dual of A/" = 1 gauge theories in four dimensions. It turns out that 
the embeddings we will find can be straightforwardly smeared in their transverse directions 
without breaking supersymmetry. Moreover, one can combine them in a way compatible with 
our generalized metric ansatz. We will use this fact to compute the backreacted geometry. 
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As the flavor branes act as a source of tlie RR forms in the backreacted solution, we will 
have to modify the ansatz of the RR three-form to include the violation of its Bianchi identity 
in a very precise form. After this modification of the ansatz, we will look again at the BPS 
equations that enforce supersymmetry and we will get a system of differential equations that 
generalizes the one found for the unfiavored system. These equations depend now both on 
Nf. and Nf and can be solved by imposing regularity conditions that are similar to the ones 
used for the unfiavored case. By solving the BPS equations for different numbers of colors 
and flavors we will discover that the system behaves differently depending on whether Nc 
is larger or smaller than ^Nf. The most interesting case occurs when A^^^ > 2A^j. In this 
regime the behavior of the solution is compatible with having an asymptotically free gauge 
theory with dynamical massless quarks. We will confirm this result by computing, from 
our solution, the beta function and the quark- ant iquark potential energy. We will get the 
expected linear confining potential and the dual description of the confining string breaking 
due to pair creation. On the contrary, when A^^ < 2Nf the solution ceases to exist beyond 
some value of the holographic coordinate. This behavior is compatible with having a Landau 
pole in the UV. 

This paper is organized as follows. In section [2] we will formulate our generalized ansatz 
for the unfiavored case. The corresponding BPS equations are obtained in appendix El It 
turns out that these equations admit a truncation, in which some functions of the ansatz are 
fixed to some particular values and, as a consequence, the system of BPS equations greatly 
simplifies. Due to this simplification we will first analyze this truncated system in subsection 
12.11 In subsection 12.21 we will consider the full system which, in general, presents a better IR 
behavior. In this subsection we carefully examine the regularity conditions to be imposed 
on the solutions of the BPS equations. 

In section |3] we consider the addition of flavor branes. We first determine the kappa sym- 
metric cylinder embeddings and then we find the particular distribution of them that pre- 
serves supersymmetry and is compatible with our metric ansatz. This distribution dictates 
the modification of the ansatz of the RR three-form needed to encompass the modification of 
the Bianchi identity induced by the flavor branes. The corresponding BPS equations for this 
case are also found in appendix |X] while, in appendix[B]we verify that, quite remarkably, the 
first-order equations derived from supersymmetry imply the fulfillment of the second order 
equations of motion for the coupled gravity plus branes theory. It turns out that the BPS 
system with flavor admits the same truncation as in the Nf = case. We study this trun- 
cated system in section |H The full system for Nc > 2Nf is analyzed in section [5l whereas 
section [6] is devoted to the study of this same system when Ac < 2A'/. 

Finally, in section [7] we recapitulate our results and discuss some possible extensions of 
our work. 
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2 Deforming the unflavored solution 



Let us begin by describing in detail the ansatz that we will adopt for the unflavored back- 
grounds we are interested in. As a particular case the family of our solutions will include 
the one found originally in [18j and interpreted in [19j as a supergravity dual of A/" = 1 super 
Yang Mills theory in 2+1 dimensions. More concretely, let cr* and a;* = 1, 2, 3) be two sets 
of SU(2) left-invariant one forms, obeying: 



da' = eijk A , dJ = ~ e^^fc A uj'' . [2.1] 



The forms a* and u' parameterize two three-spheres. In the geometries we will be dealing 
with, these spheres are fibered by a one- form A\ The corresponding ten-dimensional metric 
of the type IIB theory in the Einstein frame is given by: 



ds^ = 



dxl, + dr' + ^{ar + ^{u;'-Ar 



(2.2) 



where dxl 2 is the Minkowski metric in 2+1 dimensions, r is a radial (holographic) coordinate 
and /, g and h are functions of r. In addition, the one-form A' will be taken as: 

A' = a' . (2.3) 

with w{r) being a new function of r. The backgrounds considered here are also endowed 
with a non-trivial dilaton and an RR three- form F3, which we will take as: 

^ = ~^{uj^ - B^) A {uj^ - B^) A {uj^ - B^) + ^F' A {uo' - B') + H , (2.4) 

where B' is a new one-form and are the components of its field strength, given by: 

F' = dB" + ]^tijkB^ hB^ . (2.5) 

In (12. 4p is a three-form that is determined by imposing the Bianchi identity for F3, namely: 

c/Fg = . (2.6) 

By using (12.11) one can easily check from the explicit expression written in (12.41) that, in 
order to fulfill (12. 6p . the three-form H must satisfy the equation: 

dH = -F' AF' . (2.7) 

4 ^ ^ 

In what follows we shall adopt the following ansatz for B': 

B' = l±lM a' , (2.8) 
2 
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where 7(r) is a new function. After plugging the ansatz of written in (I2.8P into (12. 5p . one 
gets the expression of in terms of 7(r), i.e.: 

= ^drAa' + e^.j, A , (2.9) 

where the prime denotes the derivative with respect to the radial variable r. Using this result 
for in (12. 7p one can easily determine the three- form H in terms of 7. Let us parameterize 
H as: 

H = ^l^/Hir)e..,,a'Aa^ Aa' . (2.10) 

Then, by solving (12.71) for H, one can verify that 7i(r) is the following function of the radial 
variable: 

H = 27^-67 + 8k , (2.11) 

with K being an integration constant. 

In the particular case in which the function g is constant and the fibering functions w 
and 7 are equal our ansatz reduces to the one considered in refs. [HI |19]. Actually, we 
will verify that the BPS equations fix, in this case, the constant value of g to be e^^ = Nc- 
Moreover, this type of solution is naturally obtained by considering a fivebrane wrapped 
on a three-sphere in seven dimensional gauged supergravity. This three-sphere of the seven 
dimensional solution is just the one parameterized by the (T*'s, while = is the SU{2) 
gauge field of the gauged supergravity and H the corresponding three-form. The expressions 
(12. 2p and (12. 4p for the metric and RR three-form of our ansatz are just the ones that are 
obtained naturally upon uplifting the solution from seven to ten dimensions. Notice that Nc 
characterizes the flux of the RR three-form and it corresponds to the number of colors on 
the gauge theory side, whereas the constant k of (12.110 is, in the analysis of [19], related to 
the coefficient of the Chern-Simons term in the 2+1 dimensional gauge theory. 

By requiring that our background preserves some fraction of supersymmetry we arrive at 
a system of first-order BPS equations for the different functions of our ansatz. In its full 
generality this analysis is rather involved and it is presented in detail in appendix |Al Let 
us mention here that the number of supersymmetries preserved by our solutions is equal 
to two, which is the right amount of SUSY expected for an A/" = 1 gauge theory in 2+1 
dimensions. Moreover, supersymmetry imposes the following relation between the dilaton (f) 
and the function / appearing in the metric (12. 2p : 

= 4/. (2.12) 

2.1 The truncated system 

As mentioned above the equations imposed by supersymmetry on the functions of our ansatz 
are obtained in appendix \M By inspecting these equations one can check that they admit 
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solutions in which the fibering functions w and 7 vanish, as well as the integration constant 
K, namely: 

w = 7 = 0, K = 0. (2.13) 

By performing the truncation fl2.13p the first-order BPS system simplifies drastically. Ac- 
tually, one can verify that it reduces to the following three equations for = 4/, h and 
9- 



-2h 



h' 



(2.14) 



To integrate this system, let us consider first the possibility of having solutions with g 
constant. It follows from the equation for g' written in (12.141) that g must be such that: 



= 



(2.15) 



Plugging this value of g in the second equation in (I2.14p one easily shows that the equation 
for h becomes: 

h' = , (2.16) 

which can be integrated immediately as: 



2WNrr 



(2.17) 



Using these values of g and h the dilaton can be readily obtained from the first equation in 
fl2.14p . namely: 



,20 



2r 

I ' 
7*4 



(2.18) 



where 0o is a constant. The solution given by eqs. (12.150 . (12.170 and (I2.18P was obtained in 
pi] as the background generated by fivebranes wrapped on a three-cycle of a manifold of G2 
holonomy. The corresponding metric is singular at r = 0. Notice also that the dilaton (12.180 
grows linearly with the holographic coordinate for r ^ 00, as it should for a background 
created by fivebranes in the near- horizon limit. 

In order to study the system (12.141) in general and find other classes of solutions, let us 
define a new radial variable p as: 

(2.19) 
(2.20) 



,2h 



and a new function F as: 



p = e 
F = e^^ 
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We will consider F as a function of p. From the equations for h' and g' written in (12.141) we 
get the following equation for F{p): 



dF 
dp 

while the equation for the dilaton is: 



2 - 



2p 



dp F{F + iV, 



F-N, 
F + N, 



1-^ 
Ap\ 



(2.21) 



(2.22) 



From the equation for h in the system (12.14p it is straightforward to verify that the jacobian 
of the change of radial variable is: 



dr 
dp 



F + N, ' 



(2.23) 



and, thus, one can write the metric as: 



ds' 



62 



(F + iV,)2 



(dpf + ? {ar + 



UJ 



(2.24) 



By inspecting (I2.2ip we recognize our special solution (12. 151) - (12. 181) as the one that is obtained 
by taking F = Nf.in (I2.2ip and (I2.22p . Another case in which the BPS equations can be 
solved analytically is when Nc = Indeed, in this case the RR three-form vanishes, the 
dilaton is constant and (I2.2ip becomes: 



^ + - = 2 

dp 2p 



(iVc = 0) . 



(2.25) 



The general solution of (12.250 can be found easily: 

F 



4 c 
o P + — 

O p2 



(iVc = 0) 



(2.26) 



where c is an integration constant. Let us write the form of this solution in a more suitable 
form. For this purpose it is convenient to perform a new change in the radial variable, 
namely: 

P = \t\ (2.27) 



and to define the constant a, related to the integration constant c in fl^:^ as c = -4aV9y3. 
In terms of these quantities the metric (12.240 becomes: 



ds^ = 6 2* 



doc 2 ~\~ 



{dr) 



K 12 



a 



* \ 2 



(2.28) 
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where 0* is the constant value of the dilaton. Notice that (12.281) is just the metric of the 
direct product of a 2+1 Minkowski space and a manifold of G2 holonomy. This metric of 
G2 holonomy is just the well-known Bryant-Salamon metric |27], which has the topology of 
IR^ X and asymptotes to a G2 cone for large values of the radial coordinate r. Notice that 
r > a in f l2.28p and as r — a one of the two three-spheres shrinks to a point while the other 
remains finite. When a = this manifold is singular at the origin r = 0. This singularity 
is cured by switching on a non-zero value of the parameter a, in a way very similar to that 
which happens to the resolved conifold. 

Having obtained the previous solutions for near-horizon fivebranes and (resolved) G2 cones 
without branes, it is quite natural to look at solutions with RR three-form whose metric 
becomes in the UV the direct product of 2+1 Minkowski space and a G2 cone. In a sense 
these solutions would correspond to going beyond the near-horizon region of the fivebrane 
background. In terms of the variables F and p it is clear that we are looking for solutions 
such that: 

F~^p + ■■■ , (p^oo) . (2.29) 

Notice that, when (12.291) holds, 1 — SF/Ap ^ for p ^ 00 and, therefore, eq. (12.221) shows 
that the dilaton is stabilized in the UV, i.e. (p constant for large p, in contrast to what 
happens in (I2.18p . Actually, one can show that for large p the solution of the differential 
equation (12.211) that behaves as in (12.291) can be approximated as: 

F = ^p-4iV, + ^-^4 + ... , (p-00). (2.30) 

3 p 4 

By plugging this expansion in the equation (I2.22p for the dilaton, one gets the following UV 
expansion: 

d4> 27Nl 1 SlNl 1 , , 

which can be integrated as: 

27iV2 1 

= 0OO - - + •••, (P ^ 00) , (2.32) 

32 p'^ 

where 0oo is the UV value of the dilaton. For small p one gets two possible consistent 
behaviors, namely F ~ p^2 , pi . Notice that in one case F diverges at small p, while in the 
other it remains finite. Actually, when F diverges at p — ^ one can show that F{p) can be 
expanded in powers of ^p as follows: 

F = + ^ {- + -^)p+ ■■■ , (p- 0) , (2.33) 

where cq is a non-zero constant that must be taken to be positive if we want to ensure that 
F > 0. Plugging the expansion (12.330 into the right-hand side of (I2.22p . one can get the IR 
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expansion of the dilaton 



3N, 



dp 



4cov^ 



+ 



9Nl 



(P - 0) . 



(2.34) 



Notice that is regular as p — > 0, although d(f)/dp diverges. In figure [T] we have plotted the 
numerical results for F{p) and 0(p) for two different values of the constant cq. 





Figure 1: F and (j) for two different values of the constant cq of eqs. ()2.33p and fl2.34p (cq = 5 
and 10) and N^ = l. 

Let us now consider the case in which F is regular as p — > 0. Apart from the solution 
in which F = for all p, there are other solutions where F is not constant and can be 
expanded near p f» as: 

F = 6ov^-2(2 + ^)p + ^^^L_ + ... , (p^ 0). (2.35) 

Notice that in this case F vanishes at p — 0. By taking 6o > one can make F positive 
for small values of p. However, one can check by numerical integration that, after having a 
maximum the function F{p) starts to decrease and becomes negative as p increases. Due to 
this pathological behavior we will consider this solution as unphysical. 

2.2 Analysis of the general system 

Let us now come back to the general ansatz fl2.2p - fl2.lip and perform an analysis of this 
system by using the new radial variable p and the function F defined in f l2.19p and fl2.20p . 
The corresponding BPS equations are written in appendix |Xl From eqs. flA.27p and flA.28P 
it is easy to verify that the equation that determines the function F[p) is given by: 

^ = . (2.36) 

dp Dp + Dp 
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where the functions A{p), A{p), D{p) and D{p) are: 

F 1 r 



A= 2-(l-«;2)^ 
A = 2N^ {w - 7) 



2p 
F 



F - Nr 



+ NcW {w — -y) — 



D = {1- w^) (F + N^) + 2Nc w{w--f) 



D 



with V being the following function of w, 7 and k: 



V 



{w - 37) (1 - w^) - Aw + 8k . 



(2.37) 



(2.38) 



The quantities and (5 appearing in f l2.36p characterize the dependence of the Killing spinors 
on the holographic coordinate (see appendix [A]). They can be written in terms of an angle a 
as f3 = cosa, (3 = sina. Alternatively, one can write tana = A/A as in flA.25p . The explicit 
expressions of A and A are given in ( 1A.24I) . In terms of the variables p and F, they are: 



(2.39) 



Similarly, the equations that determine the functions w{p) and 7(p) can be easily obtained 
from flA.26p and flA.lSp . Let us write them as: 



dw _ B(3 + B/3 
dp ~ D(3 + D(3 



d-f _ C/3 + C/3 
dp ~ D(3 + iD(3 



(2.40) 
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where D{p) and -D(p) are the same as in (12.371) and the new functions B{p), B{p), C{p) and 
C{p) are: 



B 



B 



C 



C 



2Nr 



V_ _ g w - 7 
4p F 



2N, 



V 
Ap 



F + 3(u;-7) 



(2.41) 



Moreover, from ( ]A.27I) one easily gets that the jacobian for the change of the radial variable 
is: 

dr y/F , , 

— = . 2.42 

dp Dp + Dp 

In terms of these quantities, the metric can be written as: 
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dxl 2 + 



(Dp + Dp) 



F 



A' 



(2.43) 



Similarly, from (lA.lSP one can obtain the differential equation that determines the depen- 
dence of the dilaton (p on the variable p, namely: 



rf0 _ Ep + Ep 
dp ~ Dp + Dp ' 

where the new functions E{p) and E{p) are given by: 

1 3(l+w;2-2w7) 



(2.44) 



E = Nr 



E 



V 

8p V P 



4p 

3{w — 7) 
2y7> 



(2.45) 



As a check of eqs. fl2.36l) - fl2.45l) one can verify that they reduce to the ones of the truncated 
system when w = 7 = k = 0. Notice that in this case P = 1 and P = V = and, as 
a consequence B and C vanish and eq. (12.401) is solved by the truncated values (I2.13p . 
Moreover, one easily demonstrates that, in this case, (12.361) and (12.441) reduces to (I2.2ip and 
(I2.22p respectively. 
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2.2.1 Initial conditions 



Given a set of initial conditions for the functions F,w,'^ and 0, and a value of the integration 
constant k, the system of equations fl2.36p . (12.401) and (12.441) can be numerically integrated. 
Let us see how one can determine these initial data in a meaningful way. First of all, let 
us fix the value of the function w{p) at p = 0. Recall (see (12.31) ) that w parameterizes 
the one-form A* which, in turn, determines the mixing of the two three-spheres in the ten- 
dimensional fibered geometry. The curvature of the gauge connection A* (defined as in (12.51) 
with — * y4*) determines the non-triviality of this mixing. Indeed, if it vanishes the one- 
forms are a pure gauge connection that can be taken to vanish after a suitable gauge 
transformation. In this case one can choose a new set of three one-forms in which the two 
three-spheres are disentangled in a manifest way. On the other hand, from the wrapped brane 
origin of our solutions, one naturally expects such an un-mixing of the two S^^s to occur 
in the IR limit p = of the metric, where it should be possible to factorize the directions 
parallel and orthogonal to the brane worldvolume in a well-defined way. Moreover, by a 
direct calculation using (12. ip it is easy to verify that for w = 1 the curvature of the one-form 
vanishes and, thus, is pure gauge. Thus, it follows that the natural initial condition 
for w{p) is: 

w{p = 0) = 1 . (2.46) 

Let us now fix the value of the constant k, by adapting the procedure employed in ref. [19] 
in the case of backgrounds that are obtained by uplifting from seven-dimensional gauged 
supergravity. In this reference the authors determined n by imposing the vanishing at the 
origin of the pullback of the RR three-form on the three-cycle of the seven dimensional 
geometry which, in our notations, is the one parameterized by the one-forms a*. In the seven 
dimensional approach this three cycle shrinks at the origin and can be naturally interpreted 
as the one on which the fivebranes are wrapped. This procedure is possibly ambiguous when 
one tries to apply it in the ten-dimensional geometry, where actual D5-branes live. Moreover, 
the solutions studied here cannot be obtained, in general, by uplifting from seven dimensions. 
Therefore, it is convenient to search for a way to fix k, directly in ten dimensions. 

We start by noting that the seven dimensional cycle, parametrized by the one-forms a*, 
does not shrink in the ten dimensional geometry and, thus, it does not look strictly necessary 
that the RR three- form flux vanishes on it. Indeed, it does not shrink even in the solutions 
found in [19] . We think that the relevant cycle, which should also be the cycle on which the 
branes are wrapped, is: 

j: = {u' = a'}. (2.47) 

To understand this, let us begin by pointing out that, even if the seven dimensional gauge 
field A^ is pure gauge at the origin when the initial condition ( 12.46P holds, it is not vanishing 
there. This non-vanishing of the gauge connection is the origin of the mixing among the two 
three cycles in the ten-dimensional fibered geometry. As we are going to argue, this mixing is 
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taken into account if one considers the cycle fl2.47p 0. It is indeed easy to see that that cycle 
S is actually shrinking in the full ten-dimensional geometry if some regularity conditions are 
satisfied. Let us require that the metric function F{p) approaches a constant finite value Fq 
as p —>■ 0, namely: 

F^Fo, (p ~ 0) . (2.48) 



The induced metric on S is: 

62 



ds 



4 



P 



w 



i2 



4 



cr 



J\2 



(2.49) 



Obviously, due to the factor in brackets in fl2.49p . — * as p — ^ if eqs. (12.461) and (12.481) 
hold and the dilaton is finite at the origin. Moreover, in order to have a non-singular 
RR flux at the origin, one should require that F^ vanishes on E when p 0. We take this 
condition as a general criterium to fix the value of k. Remarkably, as can be easily verified 
from our ansatz, the pullback of F^ on S is independent of p and given by: 



4 y'^~2)^ 

Therefore, it is clear that we must fix the value of the constant k, to the value: 



= i . (2.51) 

Notice that this is exactly the value of k used in [19]. Let us see how one can reobtain this 
same value of k by requiring that the dilaton is finite at p = 0. Let Vq be the value of the 
function V defined in (12.381) at p = 0. Let us assume that Vq ^ and that w and F satisfy 
the initial conditions (12.461) and (I2.48p . Then, by inspecting (12.390 one concludes that A 
diverges at p ~ 0: 

/ — 1 

A «^ ^ K, ^ , (P ~ 0) , (2.52) 

while A remains finite at p = 0. This means that /3 ~ as p — and, therefore, the 
differential equation (I2.44p for the dilaton reduces approximately to: 

(P~0). (2.53) 

dp D 

Moreover, from (I2.45P and (I2.37p we get that the leading behavior of the coefficients E and 
D as p ^ is: 

E^-N.^^/F},, D^^Vo^, (p ~ 0) . (2.54) 



ip2 4 p 



)2 



^Alternatively, by performing a gauge transformation to A* one can get a new gauge connection = 
CT*, where ct* is a new set of left- invariant one- forms. In this new gauge the condition (|2.46p implies that 
A* vanishes at the origin and the analogue of the cycle S is just the cycle parameterized by the ct^'s with 
w' = 0. 
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Therefore, the first-order equation (12.531) for the dilaton becomes: 

dd) 1 



which, upon integration, gives rise to the divergent IR behaviour: 

~ logp + o(p) . (2.56) 
The only way to escape this conclusion is by requiring the vanishing of Vq, namely: 

Vo = . (2.57) 
But, from the expression for V in (12.381) . we get that: 

^o = 8(/t-^), (2.58) 

and, thus, the condition (I2.57P fixes again the value of the constant k to that written in eq. 
(I2.5ip . Notice that, contrary to what happens in (I2.52p . A does not diverge at p = when 
Vo = 0. Actually, A — > in this case and, therefore, the only possibility of having /5 ~ 
for p ^ 0, as is required to deduce (I2.53p . is by imposing that A vanishes faster than A as 
p — > which, in particular, implies that we must require: 

A(p = 0) = . (2.59) 

If, on the contrary, (12.590 is not satisfied, one has that /3 ~ as p — > and d(j)/dp E/ D 
1/ Pi which, again, gives rise to the undesired behavior (p ~ logp near p ^ 0. Thus, in order 
to have a regular dilaton at p = 0, we should impose the condition (12.590 . Actually, from the 
expression of A in (12.390 . as well as the initial conditions (12.460 and (12.480 . it is immediately 
possible to conclude that (I2.59P implies that the IR value of 7 should be fine-tuned to the 
value: 

7(p = 0) = 1 . (2.60) 

If (I2.60p holds, equation (I2.53P is still valid and one can check that, indeed, the dilaton 
remains finite in the IR. 

It is also interesting to look at the IR form of the metric (I2.43P when the initial conditions 
just found are satisfied. Since in this case /3 — >^ 0, only the behavior of D near p — >^ is 
relevant. One has: 

D ^ -2^oP^ , (P~0). (2.61) 

Using this result in (I2.43p . one gets that the (p, cr*) part of the metric near p ~ takes the 
form: 

% + f ("'? . (2.62) 
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Let us now change the radial variable to: 



p = . (2.63) 
The resulting metric in the (r, cr*) sector is: 

dr' + ^ {ar , (2.64) 

which is just the metric of flat four- dimensional Euclidean space. Thus, one expects that 
the metric for these solutions is regular at r = 0. We have verified this fact by explicitly 
computing the scalar curvature for our solutions and by checking that it remains finite at 
r = 0. 

2.2.2 Explicit solution 

Let us now solve the BPS equations in a series expansion around p ~ 0. For this purpose, 
let us suppose that F{p) is given by the series: 

F = Fo + Fip + F2p^ + ■■■ . (2.65) 

Then, by plugging this expansion into the BPS equations one can get the coefficients Fn for 
n > 1 in terms of Fq. The corresponding expression of Fi and F2 is: 

(Fo-iV,)(9Fo + 5iV,) 



12Fi 



(Fo - Ar,)(36F3 - 4F^N^+ 19FoN^ + 23N'^ 

^1 



p = -c;v — u -"--^ ' ----- ' ' , (2.66) 



Interestingly, one can verify that when Fq = the coefficients F„ vanish for n > 1 and 
the exact solution is F = A^c as in the background studied in [T9]. Similarly, for the initial 
conditions at p = displayed in fl2.46p and (12.601) . the functions w{p) and 7(p) can be written 
as: 

w = 1 + Wi p + W2 p"^ + ■ ■ ■ , 

7 = 1 + 7ip + 72 p' + ••• , (2.67) 
with the first two coefficients given by: 

2N^ - 3Fo 18N^ - 19 FqN^ - 16F^ + 18F^ 
= -^Ff- ' = 36^f ' 



4iV2 - 4FoiV, - F^ 
72 
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Figure 2: F and for three different values of the initial condition Fq {Fq = 2,5 and 10) for 
iVc = l. 

One can verify that for Fq = Nc one has Wn = 7n for all values of the index n. Indeed, 
in this case our generalized solution collapses to the solution studied in [1^. Moreover, the 
functions f3 and f3 behave near p ^ as: 

^ = ^ (3 - f^) + ^(/^^) ' ^ = -1 + o{p) , (2.69) 

while, for p — > the dilaton can be expanded as: 

= 0o + ^P + o(p2), (2.70) 

where 0o is an integration constant. In particular this result implies that is regular at 
p = 0, as claimed above. 

The system of BPS equations can be solved numerically with the initial conditions just 
found. From this numerical analysis we notice that, in addition to the solutions analyzed in 
[T9] (for which F = and w = 7) there are others which, for p — 00, behave as: 

F->cx), w^O, 7^7^, (p^oo), (2.71) 

where 7^, is a finite value. In figure [2] we have plotted the function F and the dilaton for 
several values of the constant Fq. These curves should be compared with the ones in figure 
[H The main differences are in the IR behavior of F{p), which is now finite at p = 0. In all 
these solutions the dilaton (p is asymptotically constant in the UV, in contrast with the ones 
of [TQ], for which the dilaton grows linearly with the holographic coordinate. In figure [3] we 
have represented the functions w and 7 for the same set of values of Fq as in figure El 

The behaviour (12.711) is easy to reproduce analytically by studying the system of the BPS 
equations. Indeed, if F, w and 7 behave as in (12.711) then one readily gets from (12.391) that. 
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Figure 3: The functions w and 7 for three different values of the initial condition Fq [Fq = 2,5 
and 10) and N^ = l. 

at leading order, A — >■ 00 and A — constant for large p and, as a consequence, j3 ^ 1 and 
/3 — > as p — >■ 00. Moreover, one can straightforwardly demonstrate that equation (12.361) 
determining F reduces to the one found in the truncated system in (12.251) . From the general 
solution written in (I2.26P we see that F ^ 4p/3 for p — > 00. Furthermore, one can verify 
that it is consistent to take the following behavior of as p ^ 00: 

w - , (2.72) 

P 

with b being a constant to be determined. Then, at leading order, one gets from (12.391) : 

4 

A ^ - p , 
3^^ ' 

A^^K-^-^Ar,7, ^A. , (p-00), (2.73) 

where K is the asymptotic value of V at p — 00. Taking into account the expression of V 
(eq. (I2.38P ). this value can be written in terms of the asymptotic value of 7 as follows: 

K = 8/t - 37, , (2.74) 

where we have momentarily considered a general value of the integration constant k. Notice 
also that [3 and [3 behave as: 

A* 9A2 1 

B ^ I ^ I — 

^ 2A2 32 p2 ' 

(p-00). (2.75) 
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Using this result one can write asymptotically the differential equation for 7(p) as: 



dp 



C 
D 



1 

J 27 



Consistency at leading order requires that 7* and K must be related by: 



7* 



9 



(2.76) 



(2.77) 



Taking into account the value of K written in f l2.74p . one gets the value of 7* in terms of 
the constant k, namely: 

Ok 

(2.78) 



7* 



2k 



Notice then that the asymptotic value of V is: 



(2.79) 



Let us now calculate the coefficient b that determines the asymptotic behavior of the function 
w. With this purpose, we notice that the functions B and B defined in f l2.4ip behave as: 



B = — 
6p - 



K + 97* 



B 



4 

7! 



(p 00) 



(2.80) 



Also taking into account that F ~ 4p/3, as well as eqs. fl2.72p and 02.750 . one gets that: 



dw 3 
'dp ~ 2^ 



K + 37, ) - 6 
which, for consistency with (12.720 . implies that: 

b = 3NrK . 



(2.81) 



(2.82) 



Finally, one can verify from fl2.44p that the dilaton (p reaches a constant value when p ^ 00. 

Taking into account that the regularity conditions in the IR fix k to be equal to 1/2, one 
gets that the actual values of 7* and b are: 



7* 



1 

3 



3iV, 



(2.83) 



a result which is confirmed by our numerical solutions. 

It is interesting to compare the background found here with the one obtained in [19]. The 
latter corresponds to D5-branes wrapped on a three-cycle of a manifold of G2 holonomy 
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in the near-horizon hmit which, as it should, has a dilaton which grows hnearly with the 
holographic coordinate in the UV. In our case the dilaton is asymptotically constant and 
the metric approaches that of a G2 cone as we move towards the UV, while in the IR region 
our solution is qualitatively similar to the one analyzed in |19]. It is thus natural to regard 
our solution as corresponding to D5-branes wrapped on a three-cycle of a G2 cone, in which 
the near horizon limit has not been taken and, thus, as we move towards the large p region 
the effect of the branes on the metric becomes asymptotically negligible and we recover the 
geometry of the G2 cone where the branes are wrapped. Notice that in reference [8] the 
authors found similar backgrounds for the case of D5-branes wrapped on a two-cycle. In this 
case the solutions asymptotically approach the conifold geometry. 

3 Addition of flavor 

Our main motivation to study a generalized ansatz of the form ( 12.2p -( !274l) was to explore the 
addition of unquenched flavors to the supergravity duals of A/" = 1 supersymmetric gauge 
theories in 2-1-1 dimensions. Indeed, we will show below that the backreacted flavored metrics 
that we will find can be represented in the form (12. 2p . i.e. their deformation with respect 
to the unfiavored ones of [I9] is of just the type studied in section O We will achieve this 
conclusion in three steps. First of all, we will study the problem in the approximation in 
which the flavor brane is considered as a probe in the unfiavored background. 

The appropriate flavor branes for our case are wrapped D5-branes that fill the Minkowski 
spacetime and are extended in the holographic direction. By using kappa symmetry [25] of 
the probe we will be able to find some simple configurations that preserve all supersymmetries 
of the background. In these configurations the D5-branes are extended along a submanifold of 
the internal space that has the topology of a cylinder and reaches the origin of the holographic 
coordinate. They can be used to add massless flavors to the gravity dual of [19]. Actually 
there is a continuous family of such supersymmetry preserving embeddings. In a second 
step we will determine how to combine these embeddings to produce a distribution of them 
that produces a backreaction on the background such that the metric is still of the form 
(12. 2p . In general the flavor branes act as sources for the RR fields and also modify the 
energy-momentum tensor. Due to the fact that we will consider a continuous distribution 
of D5-branes, these extra terms are not localized and, as we will see, their influence in the 
background can be obtained. 

The main modification of the backreacted ansatz with respect to the one studied in section 
|2] is that in the unquenched case the new RR source terms give rise to a violation of the 
Bianchi identity for F3. In a third step we will determine the appropriate modification of F3 
that gives rise to the desired violation of the Bianchi identity. Moreover, once F3 is known 
we can use it in the supersymmetry variations and obtain the BPS equations of the flavored 
backgrounds, exactly in the same way as in the unfiavored system. This analysis is performed 
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in appendix |A1 whereas the study of the different solutions of the BPS system will be carried 
out in the remaining sections of this paper. In appendix [B] we show that the backgrounds 
obtained in this way solve the second order equations of motion of the supergravity plus 
brane system. 



3.1 Supersymmetric probes 

Let us consider a D5-brane probe in some of the backgrounds studied in section [2] and 
let (/i = 0, ■ ■ ■ , 5) be a set of worldvolume coordinates. If X'^ denote ten-dimensional 
coordinates, the D5-brane embedding will be characterized by a set of functions X*^(^^). 
The induced metric on the worldvolume is: 

G(fj = d^X^'d.X'^GMN , (3.1) 

where Gmn is the ten-dimensional metric. The embeddings of the D5-brane probe that 
preserve the supersymmetry of the background are those that satisfy the kappa symmetry 
condition [28] : 

r«e = e, (3.2) 

where is a matrix that depends on the embedding of the probe and e is a Killing spinor 
of the background. Acting on spinors e such that e = ie* (as the ones of our background, see 
( ]A.6[) ) and assuming that there is no worldvolume gauge field, the matrix for a D5-brane 
probe is [251128]: 

where Gq is the determinant of the induced metric Cf} and 7^1. ..^^g is the antisymmetrized 
product of worldvolume Dirac matrices 7^. In order to define these induced matrices, let us 
denote by the coefficients that appear in the expression of the frame one-forms e— of 
the ten-dimensional metric in terms of the differentials of the coordinates, namely: 

= EfdX'' . (3.4) 

Then, the induced Dirac matrices on the worldvolume are defined as 

7^ = 9^X*^Ef Fjv , (3.5) 

where F^v are constant ten-dimensional Dirac matrices. Moreover, the puUback of the frame 
one-forms e— is given by: 

P[e^] = Efd^X^'d^'^ = Cfd^'' , (3.6) 
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where, in the last step, we have defined the coefficients = E^d^ X^. Notice that the 
induced Dirac matrices 7^ can be expressed in terms of the constant F's by means of these 
same coefficients C^, namely: 

7, = C^rM. (3.7) 

In order to obtain the particular class of D5-brane embeddings that we are interested in to 
add fiavor to the supergravity dual of 2+1 dimensional gauge theories, let us parameterize 
the forms a* and cj* in terms of the angles 6i, cpi and ipi [i = 1,2) as follows: 

= cosipidOi + sin-?/;! sin 6'i(i(y9i , = cosip2d92 + smip2sm92dLp2 , 

= — smtpidOi + cos-ipism9id(fi , cu'^ = — sin ■?/'2'^^2 + cos^/'2 sin6'2(i(y92 , 

= dipi + cos9id(pi , uj^ = dip2 + cos92dip2 . (3.8) 

Next, let us consider a D5-brane probe with worldvolume coordinates: 

e = {x\x\x',^,,r,^2) , (3.9) 
and let us embed the D5-brane in the general geometry in such a way that: 

9i = constant , ipi = constant , (i = 1,2) . (3.10) 

Let us choose the same vierbein basis as in (1A.2I) . Then, for the embedding (I3.10p the 
pullbacks of the frame one-forms are: 

p[e^'] = e^dxS (i = 0,l,2), P[e'] = eUr , 

P[e'] = P[e'] =0, P[e'] = ^ #1 , 

P[e'] = P[e'] = , P[e^^] =^[#2-^^#i]. (3.11) 

Therefore, the induced 7-matrices are: 

7^.i = e^T^i , {i = 0,1,2) , 7^ = e-^ , 

pf 1 _|_ 711 pf+9 

7^, = y ^3 - ^ e^r^) , 7v^. = , (3.12) 

and, thus, the kappa symmetry matrix of eq. (13. 3p is: 

T/t / 'yx^x^x^ibirih2 • (3.13) 
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By using the expression of the induced Dirac matrices written above, we get: 

r^e = — , ^ Tx0xix2TsTrT^e . (3.14) 

Moreover, in the type IIB theory the total ten-dimensional chirality of the spinors is fixed. 
Thus: 

r^oxix^ ri23 ri23 e = e. (3.15) 

Taking into account that ri2f = e (see eq. (1A.6P ). we conclude that: 

r^oxixs r^rg e = e . (3.16) 

Moreover, by computing the determinant of the induced geometry for these embeddings, we 
arrive at: 

-G, = . (3.17) 

From the last two equations it follows that, indeed, T^e = e, i.e. these cylinder embeddings 
preserve all supersymmetries of the background. 

3.2 Smeared configurations 

Notice that the embeddings just considered are mutually super symmetric for any value of 
the transverse angles 6i, ipi, 62 and !f2- Thus, if we have a stack of Nf flavor branes, with 
Nf 00, we can distribute them in an homogeneous way along the directions transverse to 
the embeddings (13.101) . As usual, the action for such a stack will be given by the sum of a 
DEI and a WZ term, namely: 



Nf 



S flavor = T,y2\ - (fie-^ \/-Ge + C, . (3.18) 



The smearing procedure amounts to performing the following substitution on the DBI term 
of (EUD: 



-T5V / d'^^e^ ^-Ge ^ -'^^ I a; sin sin ^2 e ^ ^ -Gg , (3.19) 



where the factor sin^i sin6'2 originates in the volume form of the space transverse to the 
embedding and (47r)^ is a normalization factor that ensures that the total number of D5- 
branes is just Nf. Similarly, the WZ term of the system of smeared flavor branes is: 

f ^^^-IT^I Vol(3^1'^) ACe, (3.20) 
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and the minus sign is due to the different orientation of the worldvolume coordinates (13.91) 
and those of the ten-dimensional space. In f l3.20p Vol( 3^]'^ ) is the volume form of the 
four-dimensional space spanned by the directions 1, 2, 1 and 2, namely: 

Vol( yl'^ ) = sin 6*1 sin 6*2 d9i A dipi A (102 A dLp2 = A A uj^ A . (3.21) 

The cyhnder embeddings just considered are extended along the 3 and 3 directions. How- 
ever, there is nothing special in our background about these directions. Indeed, both in the 
metric and in the RR three-form, we are adopting a round ansatz which does not distinguish 
among the directions of the two three-spheres. Actually, by using an appropriate coordinate 
parameterization of the a* and cj* one-forms one can straightforwardly construct supersym- 
metric cylinder embeddings that span the 1, 1 or 2, 2 directions. The volume forms of the 
spaces transverse to these embeddings are clearly: 

Vol( yl'^ ) = Aa^ Auj^ Auj\ Vol( 3^]'^ ) = a'^ A A co^ A lo^ . (3.22) 

To construct a backreacted supergravity solution with the same type of ansatz as in (12. 2p we 
should consider a brane configuration that combines these three possible types of embeddings 
in an isotropic way. The corresponding transverse volume form of this three-branch brane 
system would be: 

yoi{y^) = voi(yl'') + voiiyl'^) + voi(3^1'') . (3.23) 

Notice that the WZ term of the action of the flavor branes can be written as: 

S^Lr = n [ nAC,, (3.24) 
where Q is the following four-form: 



Mio 



^ = -TH^ Vol(3^4) . (3.25) 

Let us now write the DBI action for the D5-brane in terms of the smearing form Q. First 
of all, we notice that Q is the sum of three decomposable pieces, namely: 

n = fi(^) + + ^(3) ^ (3_26) 

where is the transverse volume form of the i^'^ branch. Let us define the modulus of any 
of these components as: 



^^....M,^....N, 11 G^^^^^. (3.27) 



k=l 
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In order to compute these moduli, it is convenient to express the f2*^'^'s in flat components 
with respect to the basis of one-forms ( ]A.2I) : 



fid) = e-4/-2/.-25 gl 



TT 



/ g-4/-2h-2g ^1 



A A A e" 



fi(3) = e-4/-2h-25 g2^g3/\g2/^ 



77^ 



It follows from the previous expressions that: 



AT 



/ „-4/-2/i~25 



(3.28) 



(3.29) 



The DBI action for the first branch in the standard coordinate system can be written in 
terms of | fi'-^-' |. Indeed, one can prove that this action is given by: 



(47r)2 



d'^'x sin^i sin ^2 6 2 \l -G^ = -T5 

Mio -'Mio 



(3.30) 



It is now clear how to generalize this result to include the three branches, namely: 



oDBI 
'-^flavor 



(3.31) 



Thus the total action of the brane distribution can be written in terms of the four-form VL. 



3.3 The backreaction 

Let us consider the coupled gravity plus branes system with action: 

S = Sub + Sfi 

avor ; 



(3.32) 



where Sub is the action, in the Einstein frame, of type IIB supergravity for the metric, 
dilaton and RR three-form F^: 



S 



1 



IIB 



2^10 



1 



1 

12 



(3.33) 



and S flavor is the action for a set of smeared flavor D5-branes, given by: 



^flavor 



-T5 [ d'\ 1^^'^ + ^5 / 



T5 / f^ACfi 
'A4io 



(3.34) 
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In (13.341) Q is the four-form that encodes the RR charge distribution of the smeared stack of 
D5-branes, while the moduh | fi^*^ | of its decomposable parts determine the mass distribution 
of the stack. In order to determine how the smeared action fl3.34p for the flavor branes affects 
the equations of motion of the RR forms, it is convenient to recall that, in the Einstein frame, 
the field strength Fy = dC^ is related to F3 as Fj = —e'^ *F3. Then, the equation of motion 
of Cq derived from the action fl3.32p is just: 

d(e-^*F7) = -2T5KjQn. (3.35) 



Using the fact that, in our conventions: 

7^5 = ^ , 2< = {2ny , (3.36) 

we can rewrite the equation for Cq as: 

d(e-'l'*Fy) = -An'^n . (3.37) 



Since F3 = —*Fy, this equation is equivalent to the following violation of the Bianchi 
identity of F3: 

dFs = in^n = -^Vol(3^4) , (3.38) 

where Vol(3^4) has been written in (I3.23p . It is clear from this last equation that, in order 
to find a solution including the backreaction of the smeared flavor branes, we must modify 
our ansatz for the RR three- form F3. Actually, we shall try to find a solution in which 

= -^3 + /s , (3.39) 

where JF3 is just given by the same ansatz (12.41) as in the unflavored case (with dj^^ = 0) and 
fs is a new term that gives rise to the violation (13.380 of the Bianchi identity. It is easy to 
see that one can take: 

/s = ^ei,fc(^* -^)^ A . (3.40) 
o Z 

By plugging the modifled ansatz (I3.39p - (l3.40p for the RR three-form into the supersymmetric 
variations of the dilatino and gravitino of type IIB supergravity one obtains a system of flrst- 
order BPS equations for the different functions of the ansatz. The corresponding calculations 
are presented in appendix |Al In appendix [B] we check that any solution of the BPS equations 
solves the equations of motion. 

Let us point out that, as happened for the unflavored case, one can consistently truncate 
the BPS equations by taking w = = k = 0. In some cases this truncation represents the 
UV limit of the solutions of the full BPS system of equations. In the next section we will 
study in detail these simplifled solutions and we will get some interesting information about 
the corresponding gauge theory duals. The analysis of the complete BPS equations will be 
performed in section [5l 
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4 The truncated system with flavor 



In this section we will analyze the truncation of the general system of BPS equations that 
corresponds to taking w = 'y = a = V = 0. In this case the equations of appendix |X] for the 
dilaton and for the remaining functions h and g of the metric reduce to: 

g' = - - N.e-^^ + ~ e"^-^^ . (4.1) 

By inspecting the system fl4.ll) one readily realizes that there are some special solutions for 
which the metric functions h and g are constant. Actually these solutions only exist when 
iVc < 2Nf and the expressions of g and h are the following: 

e^^ = 4iV, - AT, , ^'' = i^^3^' (A^. <2Ar,), (4.2) 

while the dilaton grows linearly with the holographic coordinate r, namely: 

Let us next consider solutions for which the function h is not constant. In this case we 
can use p = e^^ as radial variable as in fl2.19p and one can define the function F{p) as in 
(I2.20p . The BPS equation for F{p) is now: 

dF (F-Ar,)(2-f)-^F 



dp F + N^- ANf 

while the equation for the dilaton as a function of p can be written as: 



d(j) A'c 



dp F{F + Nc- 4A^ 



f) 



4p V 



(4.4) 



(4.5) 



Moreover, from the second equation in (14.11) we can obtain the relation between the two 
radial variables r and p, namely: 



± = VZL . (4.6) 

dp F{p) + N^-ANf ^ ' 

Notice that, contrary to the unflavored case (see eq. f l2.23p ). the sign of the right-hand 
side of (14.61) could be negative when A^/ ^ 0. This means that we have to be careful in 
identifying the UV and IR domains in terms of the new radial variable p. Let now study 
different solutions of eqs. (I4.4p - (l4.5p . 
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4.1 Linear dilaton backgrounds 



It is clear from (14.41) that in this case F = is no longer a solution of the equations. 
However, there are solutions for which this constant value of F is reached asymptotically 
when p oo. Indeed, one can check this fact by solving (14. 4p as an expansion in powers of 
1/p. One gets: 

F = N, + N,Nf^ - ^N,Nf{N,-ANf)^ + ■■■ , (p ^ oo) . (4.7) 

By plugging the expansion (14.71) into (14. 5 p one can prove that, when Nc ^ 2A^/, these 
solutions have a dilaton that depends linearly on p in the UV and, actually, one can verify 
that: 

^ - + • • ■ , (p ^ oo) . (4.8) 



dp 2{N,-2Nf) 8(7V, - 2iV;))2 p 

Notice the different large p behavior of the dilaton in the two cases N^. > 2Nf and N^. < 2Nf. 
Indeed, when A^^^ > 2Nf the dilaton grows linearly with the holographic coordinate p (the 
behavior expected for a confining theory), while for Nc < 2Nf the field decreases linearly 
with p. This seems to suggest that the beta function of the dual gauge theory depends on 
Nc and A'^^ through the combination Nc — 2Nf. Actually, one can verify that when A'c > 2Nf 
the sign of dr /dp is positive, while if Nc < 2Nf the derivative dr /dp changes its sign and r 
decreases when p increases. Indeed, by plugging the expansion (14. 7p on the right-hand side 
of (14. 6 P one gets: 

_ y/jVe ^cN] 1 

dp ~ 2{Nc - 2Nf) + 2{Nc - 2Nff p ^ ' ' ' ' IP ^ ■ (4.9) 



The first term on the right-hand side of (14.90 is clearly dominant for p — >■ oo. Its sign is 
the same as the one in the combination Nc — 2Nf, which shows that, at least in the p ^ oo 
region, the relation between the two radial variables r and p is the one described above. One 
can confirm this behavior by numerical integration of the differential equations (14.40 . (14.50 
and (14.60 . In figures H] and [5] we have plotted the result of this integration for a case in which 
Nc > 2Nf, namely Nc = 4, Nf = 1. We have integrated (14.40 by imposing the behavior (14.70 
on F{p) for large p. Once F{p) is known one can obtain 0(p) and r(p) by direct integration 
of the right-hand sides of eqs. (14. 5 p and (14.60 . We notice in figure H] that F diverges for 
p — >^ 0, while the dilaton remains finite for small p. It is easy to characterize analytically 
these behaviors. Indeed, for small p it is also possible to solve the equations (14. 4p and (14.50 
in a series expansion near p 0. For the function F{p) one has: 

f (,) . ^ , 2iN. - - i^^^^ + (1 + ) p + . . . . (4.10) 
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Figure 4: The function F and the dilaton for the case N^. = 4, Nf = 1. 
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Figure 5: Relation between the radial variables r and p for the case Nc = 4, Nf = 1. 

where cq is an integration constant which, for consistency, must be taken to be positive. In 
general, only for one particular value of Cq 7^ does one get a solution that behaves as in 
(14. 7p for p ^ 00. Similarly, the dilaton ioi p ^ behaves as: 



d(j) 
dp 



3{N, - 4Nf 



1 



9{Nc - ANf) 
4r2 



+ 



(p^O) 



(4.11) 



(Compare eqs. (14.101) and (14.111) with the ones corresponding to the unflavored solutions, 
namely (12.331) and (I2.34p ). Notice that, as in our numerical integration, eq. (14. lip implies 
that (j) is regular at p ^ when cq 7^ (although ^ is divergent). One can also easily get 
the value of the derivative dr/dp for small values of p, which is given by: 

dr _ p2 2{N^-4Nf) 



dp 



3 

P4 + 



(P ~ 0) . 



(4.12) 



From the plot of r versus p of figure [5] we notice that r grows monotonically with p in this 
Nc > 2Nf case. This means that the UV region r ^ 00 corresponds to large values of p. 
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Below we will study the beta function of the gauge theory and we will conclude that the 
theory is asymptotically free when Nc > '2Nf, while it develops a Landau pole in the UV 
when Nc < ^Nf. We can confirm this statement by looking at the result of the numerical 
integration when A^^ < 2A^/. In figure [6] we present the result of this integration for A*";, = 1 
and Nf = 1. As before, we impose the behavior (14 .7^ for large values of p. We notice that 
F{p) becomes negative at some finite value of the coordinate p, which means that the 
space ends at p = and we should consider the region p > p^ a.s the one that is physically 
sensible. Actually, in this region the dilaton (p decreases with p. A glance at the r — p relation 
displayed in figure [7] shows that r decreases with p and, actually, large values of p correspond 
to small values of r, i.e. to the IR region of the dual gauge theory. It is also clear from figure 
[7] that there is a maximal value of r, which corresponds to the minimal value p* of p. This 
fact is signaling the presence of a Landau pole in the UV of the gauge theory dual. 




Figure 6: The function F and the dilaton for the case Nc = 1, Nf = 1. 
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Figure 7: The function r(p) for the case Nc = 1, Nf = 1. This curve is obtained by 
integrating eq. (14.61) in the region in which F is non-negative. 

When Nc = 2Nf the expansion (14.81) is clearly not valid and we are in a borderline case. 
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One can prove that in this case 

| = ^p-^ + -. (A'. = 2iV,.p-.oo). (4,13) 

and, thus, grows quadratically with p when p ^ oo. Similarly, when Nc = '2Nf, the 
relation between the two radial variables r and p is: 

dr 2 2 

p+ + ... , {N, = 2Nf, p^oo) , (4.14) 



dp (AT,) I VN, 

which means that r p^. By combining the last two equations we conclude that the dilaton 
grows linearly with the r coordinate also in this Nc = 2Nf case. 

4.2 Flavored G2 cone 

Let us now consider the solution of the equations (14. 4p and (14.51) that leads to a metric which 
is asymptotically a G'2-cone with constant dilaton in the UV. It can be checked that there 
exists a solution of (14. 4p which can be expanded for large values of p as: 

4 , , 15Ar2 - 39NcNf + 24Nj 

F = -p + A{Nf- Nc) + ^ ^ + ■ • • . (4.15) 

6 p 

The corresponding expansion for d(j)/dp is: 

dd) I 27 , 1 , , 

which can be integrated as: 



4 p 32 ^ ■^>2 



- - Nc{Nc + 2Nf)- + ■■■ . (4.17) 



Notice that, when A^j = 0, the expansions (14.151) and (14.161) reduce to the ones displayed 
in eqs. (I2.30p and (I2.3ip . To find the solution in the whole range of the radial coordinate 
one can integrate numerically the system (I4.4p - (I4.5I) by imposing the asymptotic behavior 
(I4.15P to the function F{p). The results for N^. > 2Nf are similar to the ones found in 
subsection 12.11 for the unflavored system and, in particular, the solution is well-defined for 
all possible values of the coordinate p. On the contrary, when A^^^ < 2Nf the solution only 
makes sense when p is greater than some po, with pq > 0. To illustrate this fact let us 
consider a particular case with Nc < 2Nf, namely Nc = Nf. In this case the BPS system 
(I4.4p - (l4.5p can be integrated analytically. We first notice that the subleading terms in (14.150 
cancel when Nc = Nf. Actually, one can check that in this case the leading term in (14.150 
is an exact solution of the differential eq. (14. 4p . namely: 

F = ^p, (Nf = Nc) . (4.18) 
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Plugging this result into the equation (14.51) for 0, one gets: 



+ log(l-^) 



{Nf = AT,) 



(4.19) 



Nf. Notice that the dilaton in 

9iV, 



AM diverges when p = p,^ 



As a check one can verify that the expansion of (I4.19P for large values of p coincides with the 
one written in (I4.17P for Nc 
where is given by: 

Q/V 

(4.20) 



p* 



To understand the origin of this divergence it is interesting to look at the change of radial 
variables in this case. Actually, when Nc = Nf the equation (14. 6 p that determines the 
function r{p) can be integrated to give: 



r(p) 



3p H — log 



2^-3v/Ar: 

L27p + 3v^ 



{Nf = Nc) 



(4.21) 



It is clear from this expression that r(p) — oo if p ^ p*. Thus, as r must be non-negative, 
we should restrict p to the range Po < p < +oo, where po > p* is determined by the condition 
r(po) = (the actual value of po depends on the value chosen for the integration constant vq 
in Km ). 



5 The general system with flavor {Nc > '^Nf) 

Let us now consider the BPS equations of appendix |X] in full generality. In this section we 
will restrict ourselves to the case N^ > 2Nf which, according to our analysis of the truncated 
system in section HI is expected to give more sensible solutions describing an asymptotically 
free gauge theory. First of all, as in subsection 12. 2[ we are going to write these equations in 
terms of the variable p = e"^^. The equation for F = e^^ as a function of p can be written as 
in (I2.36p . where now the coefficients A, A, D and D are given by: 



A 



2-(l-»^)^ 



F - Nr 



F F 

2Nf — + NcW {w - -f) — 



A = 2Nc {w - 7) 



D = (1 - w^) {F + Nc) - ANf + 2Nc w {w - -i) 



D 




+ Nciw 




2w\/ Fp 



(5.1) 
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with V being the function of w, 7 and the constant k, written in (lA.Sp . In this case we can 
also represent /3 and (3 as in ( ]A.25I) . but now the functions A and A are the following: 



A = p + 



F H ^ (1 -i + 



4 



4 



Nr 



3F 



A 




'W\/~pF + -r- — 7) 



(5.2) 



Similarly, the equations that govern w and 7 can be written as in fl2.40p . with the coefficients 
B, B, C and C given by: 



B 



B 



C 



C 



2Nr 



4 
3 

2 
3 

4 
3 



4p ^ F 



3 - 




^F + 3{w--f) 



3 
4 



1 4iV. 2 




(5.3) 



Finally, the BPS equation for the dilaton can be represented as in (12.441) . where now the 
functions E and E are: 



E = Nr 



1 3 4Ar/ 2 o 

1 + w- 2w-y 

F 4p V . 



E 




(5.4) 



As in the unfiavored case, in order to solve the BPS system we have to fix the initial con- 
ditions of the different functions of the ansatz, as well as the constant k. To determine these 
values we follow, step by step, the procedure employed in subsection 12.21 for the unfiavored 
case. Namely, we will impose certain regularity conditions in the IR. Notice that, according 
to our analysis of the truncated system in section HI we expect that for A^^, > 2Nf this IR 
region will correspond to p ~ 0. 

First of all, let us point out that the arguments given in order to arrive at eq. (12.461) 
are still valid in this case and, therefore, we will continue to require that w{p = 0) = 1. 
Moreover, the requirement that F{p) is regular at p = is also quite natural and, thus, we 
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will also assume that (12.481) holds in this flavored case. Notice that the cycle S defined in 
fl2.47p also collapses in the IR in the present case (the induced metric on S is still given by 
(I2.49P ) and, therefore, we should impose the vanishing of the corresponding pullback of F^. 
Actually, it is immediate from (13.401) that the pullback on E of the fiavor contribution to 
the RR three-form is: 

_ 3Nf 



/sis 

and thus, the total pullback of F3 is: 

^1 



1 3Nfi 



K — 



+ 



2Nr 



cr^ A A 



Thus, we will require that n takes the value: 



K 



3N 



f 



2Nr 



(5.5) 



(5.6) 



(5.7) 



The next requirement that we will implement is the regularity of the dilaton in the IR. 
Since the reasoning that leads to the behavior (12.561) is also valid for the fiavored system, we 
conclude that we should also impose (I2.57P in the present case. Moreover, by substituting 
ti; = 1 in the expression for V in flA.5p . we get: 



Va 



K 



1 3Nf 
- H 

2 2Nr 



which vanishes precisely when k is given by the value displayed in (15. 7p . As we discussed in 

at the origin. 
0. From 



subsection 12.21 the vanishing of Vq is not sufficient to ensure the finiteness of • 
Indeed, in addition, we should require ( ]2.59p . i.e. that A is also vanishing at p 
the expression of A given in (15. 2p one discovers that this condition determines the IR value 
of 7 to be: 

, 2Ar^ 



7(P = 0) 



(5.9) 



Notice that the only freedom left by our IR regularity conditions is the value Fq of the 
function F{p) at p = 0. By changing this value of Fq we can select some particular classes of 
solutions. We are mostly interested in the backgrounds for which the dilaton grows linearly 
with the holographic coordinate in the UV and such that the function F{p) reaches a constant 
value when p — >• cx). Those backgrounds are the fiavored analogue of the ones studied in 
[TU] and can be naturally interpreted as the gravity dual of 2+1 dimensional gauge theories 
with quarks transforming in the fundamental representation of the gauge group. They will 
be obtained in subsection 15.11 by fine tuning Fq to some particular value, that depends on 
the numbers of colors and fiavors. In the UV region p ^ 00 we expect that these new 
backgrounds will coincide with the solutions of the truncated system studied in subsection 
14. H while for p — i> a significant difference between the truncated and untruncated solutions 
is expected. 
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Notice that the initial conditions fl2.46p and (15.91) ensure that the (p, a*) part of the metric 
is of the form fl2.62p . Nevertheless, in this flavored case the explicit calculation of the scalar 
curvature for the linear dilaton solutions shows that the metric is singular at the origin of 
the radial coordinate. However, as our initial conditions are such that the dilaton is finite at 
the origin, the value of the gu component of the metric is also bounded and then, according 
to the criterium of [29], the singularity is "good" and the background can be used to extract 
non-perturbative information of the dual gauge theory. 

As in subsection 14.21 we will also have backgrounds such that their metric asymptotes 
in the UV to that of a G2 cone with constant dilaton. They will be briefly discussed in 
subsection 15.21 



5.1 Asymptotic linear dilaton 

As explained above, we are interested in solutions of the BPS equations such that asymp- 
totically F is constant. Actually, by solving the BPS system in powers of 1/p, one can check 
that there are solutions in which F has the following asymptotic behavior: 

F = N, + ^ + ^ + ^ + --- , (5.10) 
p p^ p-* 

where the coefficients ai, a2 and are given by: 

ai = NfNc, 



a2 



as 



NfN, 
16 



21m 



148 Nf + 240 A^' 



(5.11) 



Notice that the first two terms in (15.101) and (15. lip coincide with the one written in (14.71) for 
the truncated system. Similarly, the functions w and 7 can be represented as: 



^1 ^ 


b2 
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fes 


w = — + 
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P' 




p3 


^1 , 


C2 


+ 


cs 


7 = — + 






P 


p2 




p3 



+ 



+ 



(5.12) 
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where the coefficients bi and Cj are the following: 



h = c, = - (iV, - 3Nf) 



b, = c2 = - (AT, - 3Nf) (AT, - 2Nj) 



3Nf 



A9N^ - 184 NcNf + 204 Nf 



49N^ - 208 NfNc + 252 



(5.13) 



By plugging the above series for F, w and 7 in the equation for 0, one can also get the 
UV behavior of the dilaton as a power series in 1/p. Actually, the ffist two terms in this 
expansion are just the ones written in fl4.8l) . For A^^ > 2iVj this means that, asymptotically, 
the dilaton grows linearly with the holographic coordinate p as: 



2(iV, - 2Ar^) 



+ o(logp) 



00 



(5.14) 
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Figure 8: On the left we plot F as a function of the holographic coordinate p for = 4, 
Nf = 1. On the right the full (dashed) line corresponds to the function ^(7) for the same 
values of and Nf. 

In order to find numerically the solution for F, w and 7 for the full range of the holographic 
coordinate one has to match the IR regularity conditions f l2.48p . f l2.46p and (15. 9p with the 
UV behavior fl5.10p - (l5.13p . As mentioned above, the only free parameter is Fq = F{p = 0). 
We have checked that such an interpolation between the p ^ and p — » 00 behaviors is 
possible by solving the BPS system with the IR initial conditions and by applying a shooting 
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technique in which Fq is varied until we obtain a solution with F{p) ^ Nc for large p. This 
only happens when Fq is fine-tuned to a very precise value. The result of this interpolation 
for Nc = 4, Nf = 1 is shown in figure [H The plot of F{p) in this figure should be compared 
with the one in figure HI which corresponds to the truncated system. The main difference 
between the results for F{p) in these two figures is that F{p) diverges for p — > in the 
truncated system, while it remains finite in the complete solution, whereas for large p both 
solutions nearly coincide. Notice that w and 7 evolve smoothly from their initial values at 
p = to their vanishing asymptotic values for large p. The dilaton, which is not shown in 
figure [8], grows monotonically with p and becomes approximately a linear function of the 
holographic coordinate when p is not very small. These features are in agreement with the 
expectation that these solutions of the complete system would reduce to the equivalent ones 
of the truncated ansatz in the UV. Notice that in the borderline case = 2Nf the initial 
value of 7 is zero. The result of our numerical calculation shows that, in this case, the 
function 7 becomes negative for p > and approaches its asymptotic vanishing value for 
p 00 from negative values, in agreement with expansion written in fl5.12p . 

Having obtained this solution of the equations of motion of the gravity plus brane system, 
let us see if it incorporates some of the features that the supergravity dual of 2+1 dimensional 
gauge theory plus fiavors should exhibit. First of all, in the next subsection we will give a 
prescription to evaluate the gauge coupling and we will verify that, for A^^^ > 2Nf, this 
coupling displays the expected property of asymptotic freedom in the UV. Moreover, in 
subsection 15.1.21 we will analyze the potential energy for an external quark-antiquark pair 
and we will discover that this potential behaves in the way expected in a theory which has 
string breaking due to pair production of dynamical massless quarks. 

5.1.1 The Yang-Mills coupling and the beta function 

Let us study the evolution of the gauge coupling constant with the holographic coordinate. 
In order to do that, let us consider a D5-brane probe extended along the three Minkowski 
directions and wrapping some internal three-cycle at a fixed value of the holographic coordi- 
nate. The natural three-cycle to compute the Yang- Mills coupling is just the one used above 
to fix the constant k, namely S = {cu* = a*}. Indeed, as shown in section |2l S shrinks to zero 
size at p — >■ 0, which corresponds to the IR of the gauge theory where one expects to have 
a 2+1 dimensional behavior of the D5-brane probe. Thus, S is the analogue in the present 
case of the two-cycle found in ref. [30j for the background dual to A/" = 1 super Yang-Mills 
in four dimensions. 

The DBI action for such a probe in the Einstein frame is: 




(5.15) 



where Gq is the induced metric on the D5-brane worldvolume and J-' is the worldvolume gauge 
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Figure 9: The inverse of the Yang-Mills coupling as a function of the holographic coordinate 

p for the case N^. = 4, Nf = 1. 

field. By looking at the JF^ terms in the above action, we get the value of the Yang-Mills 
coupling constant of the dual 2+1 gauge theory, namely: 



1 -3, 
= p 4 ' 



9ym 



^ y^-det(G3) d'^, (5.16) 



where the induced metric on the three-cycle E has been written in eq. fl2.49p and we have 
neglected all constant numerical factors. By using the metric written in (12.491) . we obtain: 



1 



9ym 



(5.17) 



where again we have neglected all numerical multiplicative constants. Due to our boundary 
condition (12.461) . the right-hand side of (I5.17P vanishes for p = 0, which corresponds to 
having — oo in the IR, as expected in a confining theory. Clearly, l/gyM grows as we 
move towards the UV region p — > oo, in agreement with the expected property of asymptotic 
freedom. In figure [9] we have plotted I/qym fo^^ = 4 and Nf = 1. In order to obtain the 
corresponding beta function from (I5.17P we would need the relation between the coordinate 
p and the energy scale of the problem. The usual arguments employed for the gravity duals 
of four-dimensional Af = 1 gauge theory are not valid in our three-dimensional case and, as 
a consequence, such an energy-radius relation is lacking here. The best that we can do is to 
use the original radial variable r which, as argued in subsection 14.11 grows in all cases when 
we move towards the UV. This fact can be further verified by placing a fundamental string 
stretched along the radial direction and looking at its energy as given by the Nambu-Goto 
action. Clearly this energy grows in the direction of increasing dilaton, i.e. when r increases. 
When Nc > '^Nf, large p corresponds to large r and, as F ^ and w ^ for our solutions, 
one approximately has: 

d f 1 \ 3 1 dp 

^ - p2 — , (p^ oo) . (5.18) 



dr V J 2 dr 
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Moreover, when > 2Nf we get from the asymptotic expansion (14. 9p : 

(p^oo), (5.19) 



dp 2{N^ - 2Nf) 



3 



dr ^/Nl 
and, thus, we can write for large r and p: 

df ^ \ {N^-2Nf) 1 {N,-2Nf)^ i 

Eq. (15.201) shows that I/qym grows with r in the UV (as it is obvious from figure [9]) and 
that its derivative also grows with r as ^/r. This behavior is consistent with the expected 
negative beta function for gyM- The form of eq. (15.201) could lead to the conclusion that its 
right-hand side vanishes in the borderline case Nc = 2Nf. However, one should be careful in 
this case and use the correct expression (I4.14p in (I5.18p . One gets for large p: 

(iVc = 2iV^), (5.21) 

which shows that, actually, that the right-hand side only vanishes when r — > oo and the 
theory is still asymptotically free. 



5.1.2 Wilson loops 

In order to verify how the flavor degrees of freedom are encoded in our backreacted geometry, 
let us study the rectangular Wilson loops for external, non-dynamical, heavy quarks. These 
Wilson loops can be evaluated by studying the Nambu-Goto action of a fundamental string 
whose ends lie in the UV region r —>■ oo and are separated by a distance L in the gauge 
theory directions [31], [32] (see also [331 El] ) • To describe such configurations let us choose the 
time t and a Minkowski coordinate x as worldvolume coordinates and let us parameterize 
the string worldsheet by means of a function r = r(x), where r is the holographic coordinate 
of (12. 2p . The induced metric in the string frame is: 

e'^l-df + {1 + {r'f)dx^] , (5.22) 

and, as a consequence, the Nambu-Goto action takes the form: 



S 



j dtdx e'f'^'^ Vl + (rO^ , (5.23) 

where we are taking the string tension to be equal to one and r ' = dr/dx. From the invariance 
of the lagrangian in (I5.23P under shifts in the coordinate x, we immediately obtain a first 
integral of the equations of motion of the string, namely: 



(r')2 



39 



where r^m is the minimal value of the holographic coordinate reached by the string world- 
sheet. From fl5.24p we can straightforwardly obtain r', with the result: 



r' = ±- ^ . (5.25) 




Figure 10: In this figure we plot the energy E of the Wilson loop versus the length L. The 
curve on the left corresponds to the case = 4, Nf = 1. This curve clearly shows that 
there is a maximal L and that E is a double-valued function of L. On the right we plot 
the same quantity for the unfiavored case with = 4, where string breaking does not take 
place. 



It is now trivial to obtain the length L, i.e. the separation between the quark and the 
antiquark, function of the minimal value of r: 

Lir^in) = 2 / dr , (5.26) 

where r^ax is a cutoff related to the mass of the external quarks, that can be taken to be 
very large. Moreover, after subtracting the masses of the non-dynamical quarks, the energy 
of the string configuration becomes: 



r-max g20(r) fr 



max 



Ehmin) = 2 / = dr - / e2<^(^) dr , (5.27) 

which can be identified with the potential energy of a quark- antiquark pair separated by a 
distance L{rmin)- Notice that both E and L depend parametrically on r^m. By varying r^m 
we can obtain the corresponding values of E and L and extract the dependence of the energy 
on the distance. In a theory with dynamical quarks one expects the strings to elongate until 
their tension equals the mass of the lightest meson and then create a quark-antiquark pair 
and break. In an E versus L plot this behavior would correspond to having a maximal value 
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of L and a double- valued function E{L). The corresponding result for our solution when 
Nc = 4: and Nj = 1 are displayed in figure fTOl We see that the expected behavior of E{L) is 
reproduced, in a way similar to one found in ref. [8] for the background dual to an A/" = 1 
SQCD-like theory in 3+1 dimensions. To allow for a comparison with the unflavored theory 
we have also plotted in figure fTOl the result of the E{L) curve for Nc = 4: and Nf = 0. In this 
unflavored case there is no maximal value of L and, for large quark-antiquark separation, 
the energy grows linearly with L, as it should for a confining theory without screening due 
to pair creation. 

It is also possible to understand the different behaviors of the E{L) curves by analyzing 
the function that is integrated on the right-hand side of (15.261) . Indeed, one can check that 
when r and r^m are both small, the square root on the denominator in (15.261) behaves as: 

^e^'ACO - e2'^('-™-) ~ r" , (r,r„i„) ^ , (5.28) 

where a is some constant. In the unflavored case, by combining eq. (12.701) and the fact 
that p ^ r"^ for small r, one concludes that a = 1, which in turn implies that the integral 
giving L{rmin) is divergent when Tmin —>■ 0, in agreement with our numerical results. On 
the contrary, when flavors are added one can verify numerically that a < 1, which means 
that the integral (I5.26P is now convergent when r^m —>■ and there is a maximal value of L. 
This different behavior of the dilaton in these two cases is correlated with the fact that the 
flavored metric develops a (good) curvature singularity at the origin, while the unflavored 
solution is regular. It is worth pointing out that similar results have been found in . 

5.2 Asymptotic G2 cones 

When Eq takes values in a certain range, the solutions of the BPS equations reached at the 
UV have a metric which is the direct product of 2+1 dimensional Minkowski space and a G2 
cone. The solutions in this case are very similar to the ones discussed in subsections 12 . 2 . 2| and 
I4.2l and we will not discuss them further here. Let us only mention that the asymptotic values 
of E, w and 7 for p 00 can be determined analytically, in a way completely analogous to 
the one employed in subsection 12.2.21 (one has to take into account the different value (15.71) 
of the constant k in the present flavored case). One gets the following asymptotic behavior: 

E ^^p + 4iNj-Nc) + ■■■ , 

3{Nc - 3Nf) 

w ^ + ■ ■ • , p ^ 00 , 

2p 

1 Nf , , 

7«3-^ + --, (5.29) 

a result that is confirmed by our numerical calculation. 
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6 The general system with flavor (Nc < '^Nf) 



In this section we briefly describe the solutions of the complete BPS system for < 2Nf that 
have an asymptotic linear dilaton in the UV. Even if the interpretation of these solutions is 
less clear than in the Nc > 2Nf case, it follows from the analysis performed in subsection 14.11 
for the truncated solutions that the right holographic variable is now the original coordinate 
r, instead of the variable p used so far. Moreover, we learned in that subsection that, in this 
Nc < 2Nf case, the small p region should be interpreted as the UV of the gauge theory (with 
large r) and vice versa, the IR would correspond to large p and small r. Somehow when 
going from Nc > 2Nf to A^^^ < 2Nf the UV and IR are exchanged. Thus, when Nc < 2Nf, it 
is natural to search for solutions that approach the unfiavored one at the IR@ where, in terms 
of the variable p, they can be represented by the series ( IS.lOp and ( 15.12^ . Notice that when 
one uses the variable r as the independent variable one should also determine the function 
p(r) or, equivalently, h = h{r), where h{r) is the function squashing the ex' sphere in our 
ansatz (12. 2p . The differential equation determining p(r) is just the one written in (I2.42p . 

F 

1 

. 8 
. 6 
. 4 
0.2 




Figure 11: The function F{p) for Nc = Nf = 1. To obtain this curve we have imposed the 
behavior (15.101) for small r (or large p). 

The numerical integration of the BPS system for A^^ < 2A^/ with the IR initial conditions 
given by (15.101) and (15.121) shows that the function F decreases from its initial value F ^ Nc 
at r ~ until it vanishes at some finite value r* of the coordinate r. For larger values of r 
the function F becomes negative and the solution does not make sense any more. In figure 
[11] we plot F as a function of r for Nc = Nf = 1. Notice that F{r) drops very fast to zero 
as r approaches its final value r^,. By choosing appropriately the initial conditions in the 
integration of p(r) (or, equivalently, of h{r)) we can make that r = r^, corresponds to p = 0. 
Moreover, the dilaton (not shown in figure [TTj) grows linearly with r, as expected for this 
type of solutions. In these calculations we have taken the same value (15.71) for the constant 
K. Notice that, with our choice of initial conditions, the cycle S collapses at r = r*. 

^This is actually what happens in other backreacted solutions with a Landau pole, such as the ones in 
[15] for the conifold. 
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As argued in subsection 14.11 for the truncated system, we think that the proper interpre- 
tation of the point r = r* is the location of a Landau pole. To confirm this interpretation 
we have calculated the Yang-Mills coupling, as given by (15.171) . as a function of the radial 
variable r. The numerical results confirm that 1/ Qym is a monotonically decreasing function 
of r and that gyj^j — cx) as we approach the point r = r*. 

7 Summary and discussion 

In this paper we have found backgrounds which encode the effect of adding a large number 
of unquenched flavors to the gravity dual of A/" = 1 gauge theories in 2+1 dimensions. By 
using kappa symmetry, we first determined the appropriate embeddings of the flavor branes 
that preserve the supersymmetries of the unfiavored background and, then we found the 
modification of the ansatz of the RR field needed to solve the Bianchi identity of the coupled 
gravity plus branes system. We have subsequently obtained a system of first-order BPS 
equations, which we have solved with different boundary conditions. The most interesting 
solutions are those that contain a linear dilaton in the UV. When N^. > 2Nf we have argued 
that these solutions display the expected properties of a gravity dual of an asymptotically 
free theory with dynamical quarks. We have checked this fact by computing, from our 
background, the Yang-Mills coupling constant, as well as the expectation value of the Wilson 
loop. In this latter case we have explicitly verified the expected string breaking due to quark- 
antiquark pair production. We also found solutions for Nc < 2Nf and we have shown that 
they are consistent with having a Landau pole in the UV of the gauge theory. 

Let us comment on some points that, in our opinion, would need some further clarification. 
First of all, it would be desirable to have a more precise characterization of the field theory 
dual to the background found here. One could argue as in [8] and try to determine the 
IR field theory that is obtained by integrating out the massive KK fields. Following the 
reasoning of [8] we conclude that extra couplings of the fundamental matter fields (quartic 
or with higher powers) are generated. For this reason some of our results are difficult to 
check on the field theory side. 

One of the problems that would be interesting to understand from the field theory point 
of view is the dependence of the beta function on and Nf. Notice that, due to the 
low amount of supersymmetry preserved by our solution, one cannot rely on the power of 
holomorphy, which has been so useful to extract the non-perturbative structure of A/" = 1 
gauge theories in four dimensions. In a certain sense these theories are a good arena to test 
the power of holography as a tool to explore the strong coupling regime of gauge theories. 
As a first step in this direction, let us try to determine the coefficient k of the Chern-Simons 
term for our solution. In general, to obtain such a result one should be able to find the 
shift of the level due to the integration out of the KK fermions. In the presence of massless 
flavors this calculation is even more complicated because an explicit computation of the 
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Witten index is, at least to our knowledge, lacking. However, the form of the supergravity 
solution that we find seems to suggest that the Chern-Simons level k of the low-energy three 
dimensional theory, after having integrated out the KK fermions, is: 



k^^^^. (7,1) 



This would imply the Witten index for such theory would be: 



nlNr 



" c- 



with n = k — {Nc — 3Nf)/2 > 0, and being otherwise. It would be very interesting to prove 
or disprove these results from a direct calculation in the field theory. 

Let us finish this section by mentioning some further topics that one could address from 
the supergravity side. First of all, one could try to generalize our background to the case in 
which the flavors are massive. For the case of the conifold this generalization was achieved in 
[T5] by a simple modification of the RR form whose Bianchi identity is violated. One could 
try to apply a similar procedure for the setup studied here. The analysis of the meson spectra 
for the backreacted geometry is clearly another interesting problem to look at. To perform 
this analysis one can add a probe and consider its fluctuations. Presumably one would find 
the same type of problems related to the normalizability of the fluctuation modes as in other 
backgrounds generated by D5-branes and a careful treatment would be needed to extract 
the meson masses. The construction of the black hole version of our background is clearly 
of great interest since it would allow us to explore the thermodynamic and hydrodynamic 
properties of the field theory dual at finite temperature. 

For the metrics with G2 holonomy one can perform the so-called flop transformation, 
in which the two three-spheres are exchanged. An interesting problem to study is to what 
extent this transformation can also be performed in our solutions and what are the effects on 
the field theory side (see p5] for a similar study in the case of backgrounds of G2 holonomy 
without fluxes). According to the ideas put forward in |8], one expects that this discrete 
transformation is a kind of Seiberg duality, that somehow would exchange the rank and the 
level of the corresponding field theory. Finally, one could try to see if our family of unflavored 
backgrounds can be used to describe the physics of domain walls in A/" = 1 gauge theories 
in four dimensions. 
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A Appendix: Derivation of the BPS equations 

The supersymmetry transformations for the type IIB dilatino A and gravitino in Einstein 
frame, when the RR three-form is nonzero, are: 



S^, = D,e + ^ F(^),^,3 ( P^^^^'^^ - 9 S^^^ P^^^^ ) 6* . (A.l^ 



We want to solve the conditions dX = 5ip^ = for a metric given by our ansatz (I2.2p . In 
what follows we shall choose the following vierbein basis: 



f = dr 



J+a 



2 ' 



(A.2) 



where is parameterized in terms of the function w{r) as in (12.31) . The spin connection in 
the basis (lA.2p is given by: 



UJ 



UJ 



UJ 



UJ' 



UJ' 



UJ 



UJ 



V 



e-ff'e^' 



e-f (/' + h') e' - ^ e' 



w 



e"^ (/' + g') - -J e' 

Jjk 



Jjk 



1-w' 



(A.3) 



We will take the RR three-form F3 that corresponds to the general system with flavor of 
the main text, i.e. we will take F3 as given by (13.391) . The different components of this field 
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strength in the basis (lA.2p are: 



rii 

i23 

ijk 

tjk 
p{3) 



where is a function of r defined as: 



2 " '^V- iV, 

3/-3/1 



(A.4) 



W + 



(A.5) 



= -37) (1 - w^) - 4(^1 

with K being the same constant as in fl2.1ip . 

The Kilhng spinors of the background are those e for which the right-hand side of eq. 
( lA.ll) vanishes. In order to satisfy the equations 5\ = = we will have to impose 
certain projection conditions on e. These conditions are: 



le = e . 



Let us now define the following matrix 



rl23 • 



(A.6) 
(A.7) 



From the vanishing of the dilatino variation under SUSY, we get: 



4A^. 



/ 



+ w'^ — 2w'y ) e 



~g~2h 



' 2 



TiiT^e = . 



(Ai 



Let us now consider the gravitino variation. From the condition Sip^i = we obtain that the 
metric function / must be related to the dilaton as: 



(A,9) 



Moreover, from the equation Sipi = 0, after using eqs. (1A.8I) and (lA.9p . one arrives at: 



h'e 



2 



[1 - vj^) + A'e 1 



^ + - 2w-f ) e ^ 



w 



rue + 



+ 



2 L 4 



FiiF.e = 



(A.IO) 
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Similarly, the condition Si/j^ = leads to: 



oa-h 



+ 



w 



r^i e + N,{w--f) e-''^-^ r^i r, e = . 

In order to solve the above equations, we shall impose the additional projection: 

r,e= (/5 + /3rii)e , (A.12) 

where {3 and (3 are functions of the radial variable to be determined. As (F^.)^ = 1 and 
{r*,r^j} = 0, by consistency, these quantities must satisfy the condition: 



(A.ll) 



(3^ + (3^ = 1, 

and, therefore, they can be represented in terms of a single angle a as: 

(3 = cos a , P = sin a . 

Notice that the projection (]A.12p is equivalent to 

and can be solved as: 

e = e 2 li eg , 

where a and the spinor eo depend on r and the latter satisfies the projection: 

r* Co = cq. 



(A. 13) 

(A. 14) 

(A.15) 
(A.16) 

(A.17) 



We can now write the set of BPS equations corresponding to this ansatz. Let us substitute 
(1A.12I) and ( 1A.15I) in the dilatino variation (1A.8I) . By separating the terms containing the 
unit matrix from those with F^^j, and using (j) = 4/, we arrive at the following equations: 



f = ^ 
J 4 



7 



-39 



-2g+h 



- 1 



+ w"^ — 2w7 ) e 



- 1 



^ + - 2w7 je"'' 



+ 3(w - -f)e-^^-^ 

4 \ ' y 



/3 



-I 



/3. 



(A.18) 
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In order to determine (5 and /?, let us plug (1A.12P and (lA.lSp in (lA.lOp and consider the 
terms containing r^|. One gets: 



L 2 



+ 2 



[w 



Similarly, from (lA.lip one arrives at: 

1 - w^) e^-'^^ - 4e-f + - N^(l - ^ + w^- 2wj^ e-^-^^ 



(A.19) 



/3. 

(A.20) 



By substituting the expression of 7' taken from flA.lSP into flA.lQp . one gets the following 
expression of w'\ 



N, 



^Ve-^^ - Awe 



/3 - 



(A.21: 



whereas, by performing a similar manipulation to flA.20p . one can prove that: 

V 



L6^ 



-3h 



+ 



3 



(A.22) 



By eliminating w' in (IA.2ip and (IA.22p . one can demonstrate that (3 and (3 satisfy a relation 
of the form: 

A/3 - A/3 = , (A.23) 
where the functions A and A are given by: 



A = e^^ ^ ^ ^ 1 

4 4 



4A^ 



/ 



Nr 



+ w'^ — 2w7 ) — — - e 



o2h-2g 



A 



N. 



24 



(A.24) 



48 



One can solve (1A.23P for j3 and /3 as /5 oc A, /5 oc A, where the common proportionahty 
function is determined by imposing the condition + = sin^ a + cos^ a = 1 (see eq. 
f lAlSl) ). One gets: 



A 



VA2 + A2 

Let us now rewrite (1A.22I) as: 

'V 



P 



A 



VA2 + A2 



(A.25) 



= AT, 



6 



^-''^ - 2{w - 7)e- 
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(3 + 



+ 



p. 



(A.26) 



Up to now we have determined the first-order BPS equations satisfied by /, 7 and w (eqs. 
(lA.lSp and (1A.26P ). Let us now determine the equations for the remaining functions of our 
ansatz, namely h and g. By using (1A.12I) and (lA.lSp in (1A.10I) and considering the terms 
containing the unit matrix , we get: 



h' 



'1 — w )e 



P + 



1 

+ 2 



p. 



(A.27) 



Similarly, from flA.lip we obtain the following equation for h: 



^ =4 



4e-9 _ (1 - u;2)g9-2h + _ 4^ + ^2 _ 2^;^ j e^^-''^ - m,e-^' 

+N,{w --f)e-'^^-^ p. 



P + 
(A.28) 



To complete our analysis of the Killing spinor equations we should check that the variation 
of the radial component of the gravitino vanishes. Actually, one can check that this condition 
holds if the following two equations are satisfied: 

(A.29) 
(A.30) 

where a is related to P and (3 as in flA.14p . By using flA.25p one can check, nontrivially, 
that eq. flA.29p is a consequence of the other first-order equations of the system. Moreover, 



a 



-eo , 
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(1A.30P determines the dependence of eo on the radial coordinate. Indeed, it can be integrated 
as: 

eo = etr] , (A.31) 

where rj is a constant spinor. By combining eqs. ( 1A.31I) and ( 1A.16I) . one can get the exphcit 
expression of the Kilhng spinors, namely: 

e = e"^^" es 77 , (A.32) 

with r] being a constant spinor satisfying the four commuting projections: 

Tii V = '^22V = Tss V, 
it]* = T] , 

T,ri = ri. (A.33) 

The projections flA.33P show that our background is 1/16 supersymmetric, i.e. it preserves 
two supersymmetries. Notice that this is the amount of supersymmetry expected for an 
A/" = 1 theory in 2+1 dimensions. 

To finish this appendix, let us point out that the BPS equations just found are consistent 
with the truncation w = 7 = k = (which implies that V = 0). It is clear from ( 1A.29I) 
that in this case the phase a can be taken to vanish or, equivalently, /3 = 1 and (3 = 0. The 
corresponding truncated equations are much simpler than the full BPS system and will be 
studied separately in the main text. 



B Appendix: Equations of motion 



The equation of motion for the dilaton derived from the action (13.321) is: 

7^*' (G"^' V^a. - ie*F| ^ - W . (B.l) 

In order to make this equation more explicit, let us compute the determinant of the metric 
for our ansatz. The easiest way to do that is by realizing that y—G can be obtained by 
computing the wedge product of all the one-forms of the frame basis: 

e^" A • • ■ A = dx^ A ■ ■ ■ dx^ . (B.2) 

In our case, we get: 

v^= ^e^^^^^'^^^) V^, (B.3) 

where gi and ^2 are angular factors that depend on the parametrization of the two sets of 
left-invariant SU{2) one-forms. Actually, let us suppose that we represent the a"s (a;"s) in 
terms of the three angles ^1, (pi and ipi {62, 02 and tp2)- Then: 

cr^ A cr^ A cr^ = dOi A A d^i , a;^ A a;^ A a;^ = a/^ d92 A #2 A d%l)2 . (B.4) 
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For the standard election of the angular variables, one can verify that: 

J^i = sinOi , (i = 1,2) . 



(B.5) 



Using these results one can verify that the equation of motion for the dilaton (1B.1|) reduces 
to: 



30 



6 ^ .^flavor 



(B.6) 



We shall evaluate the first two terms on the left-hand side of the above equation by using 
the BPS equations for the functions of the ansatz and the dilaton. Notice that, in order to 
compute the second derivative of the dilaton, we must differentiate its BPS equation. After 
this differentiation we need to evaluate the derivatives of (3 and (3. Taking into account 
( 1A.14I1 . one gets: 

p' = -pa' , (3' = (3a' , (B.7) 

and, after using the value of a' from (lA.29p . one obtains: 



3 
4 



N.^'e-^-'' + -w'e^-'' 



(3 



p' = ~ N^ie-<^-^ + ^w'e<^-^ 



(B.8) 



Moreover, from our ansatz the term containing the RR three-form in the equation of 
motion of the dilaton is: 



3(j> 

12 ^ 



3 



1 - 



AN 



N, 



f 1 
+ to — 21^7 



+ - {w - 7)^ e" 



4g-2h 



+ 



Ar2 



2 ^ ' " 32 
Amazingly, after a very long calculation one gets the simple result: 



h" + 2(0' + 3(/.' + g')^' - ^ = 6Ar^e-2^-2^ . 



(B.9) 



(B.IO) 



In order to verify that the right-hand side of (1B.6P reproduces exactly this result, let us use 
the action of the smeared flavor branes in terms of the charge distribution four-form f2. From 
the expression of this action it is straightforward to evaluate its contribution to the equation 
of motion of the dilaton: 



2K^g 



I 7^ ^ ^ 1 ^flavor 



:b.iii 



which proves that, indeed, the equation of motion of the dilaton is satisfied. 
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Let us next verify that the Einstein equations are satisfied. These equations are: 

+ 



Rmn — -^GmnR 



^ (dM(l>dN 



^GMNdp(j)d^(, 



1 



MN 



(B.12) 



where Tmn is the DBI contribution to the energy-momentum tensor, namely: 

Tmn 



2K?f) SSdbi 



(B.13) 



Taking into account that the absolute value of the smearing form VL is defined in fl3.27p . 
we can simply compute the DBI contribution to the energy-momentum tensor, in fiat com- 
ponents: 



■ MN 



3! I m 



(B.14) 



In components, it means that 



T 

■J- rr 



2Nf e-2/-2h-29 



(B.15) 



with the off-diagonal components being zero. 

Moreover, from the value of the spin connection we can easily compute the different 
components of the Ricci tensor. In flat indices these components are the following: 



Rfi 



-Vx.x'^e-'f [f" + 8(/0' + ?>fh' + ?>rg'] 



3f" + h" + g" + h'{f' + h') + g'if' + g') + ^^^'^w')' 



/" + h" + Sif'Y + Sih'Y + llf'h' + 3/V + 3/iV + 



+ e 



2g-4h 



[l-W 



2\2 



+ e 



2g-2h 



[W 



i\2 



8 



- 2e 



-2h 



-2f+g-h 



+ 2nj'f' + -w'g' + -w'h' + e'^^-^^ ""^^ 



— e 



-2/ 



id" ^ 1 

— + 2w'f' + -w'g' + -L , _ 
4 4 ^ 4 2 

f" + / + Kff + Kg'? + 11/y + ^fh' + ?>h'g' - 



— e 



[l-W 



2\2 



=29 



8 



i\2 



2e 



-29 



(B.16) 
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with the remaining components zero. It is straightforward to calculate now the curvature 
scalar in Einstein frame: 



R 



— e 



-2/ 



18/" + 6h" + 6g" + 72{f'Y + 12(/i')' + 12(^') + 5Af'h' + 5Af'g' + 18h'g 



(B 



Finally, one can verify that the BPS equations we have found imply the fulfillment of eq. 
(IB. 121) along its different components. 

Finally, it remains to verify the equation of motion for the RR three-form F3, which is 
just: 

d{e^*F^) = . (B.18) 

After computing the Hodge dual of with the metric (12. 2p . one can demonstrate that 
( IB. 181) is equivalent to the following second-order differential equation for the functions of 
our ansatz: 

e^^+s [/ + (2c/,' + h' + g')i] + 4e'^-^(u;-7) + (l - w^)V + 



+ 2e^-^{^l - ^ + - 2^7) w = 



:B.19) 



which, again, can be shown to be a consequence of our first-order BPS equations. 
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